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Abstract 



The charm contribution to structure functions of diffractive deep inelastic 

scattering is considered here within the context of the Ingelman-Schlein model. 

Numerical estimations of this contribution are made from parametrizations 

of the HERA data. The influence of the Pomeron flux factor is analized as 

well as the effect of the shape of the initial parton distribution employed in 

the calculations. The obtained results indicate that the charm contribution 

to diffractive deep inelastic processes might be large enough to be measured 

in the HERA experiments. 
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I. INTRODUCTION 



The HERA data of deep inelastic scattering (DIS) measured in the last few years contain 
a sizeable fraction of events with a large rapidity gap in the forward region This phe- 
nomenon is present even at high and results from a colour-singlet exchange between the 
dissociated virtual photon and the recoiling proton (or proton remnant), characterizing what 
is usually called diffractive DIS (DDIS). The measurement of DDIS at HERA provides an 
unique opportunity to study diffraction in regions in which perturbative QCD is applicable. 

Open heavy quark production at HERA is also a subject of major interest in QCD 
phenomenology. Both the HI and ZEUS collaborations have found the charm component 
of the structure function, F2'^\x,Q'^), to be a large fraction of F2{x,Q'^) at small x 
Recently the first measurements of the 66-cross section have been reported 0. Due to the 
higher mass of the 6-quark, it is two order of magnitude smaller than the cc-cross section . 

For the moment most of the experimental data is for the neutral current structure func- 
tion F2{x,Q'^). A number of theoretical estimates of F2'^\x,Q'^) has recently been ob- 
tained (see also a review in Ref. and references therein). In the present paper 
we consider DDIS with open charm production and calculate a charm contribution to the 
diffractive structure function F-f {(3 , Q"^ , x-p) . 

There are two different approaches to a treatment of the charm component in structure 
functions. In one approach [^|lO|l the charm is an active flavor which undergoes massless 
renormalization group (RG) evolution. We will follow another approach in which only light 
{u, d, s) quarks and gluons are active partons and no initial state heavy quark lines show 



up in any diagram [|r^. It involves the calculation of the photon-gluon fusion process and 
thus is quite sensitive to the gluon distribution. 

For the time being, the gluon distributions inside the Pomeron predicted by a number of 
models are dramatically different and they have different shapes. The diffractive production 
of open charm at HERA will, therefore, provide us the possibility for a direct test of the 
models. 
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The present analysis is partially based on a previous study on the Pomeron structure 



function in which charm contribution was not considered. This study was mostly con- 
cerned with effects of the Pomeron flux factor on the evaluation of the diffractive structure 
fuction. Such effects are a central issue also in the present analysis. 

This paper is organized as follows. In Section II, we describe how to take into account the 
charm content of the Pomeron. The charm contribution to the diffractive structure function 
is evaluated in Section III, where the comparison with other models is also given. Our main 
conclusions are summarized in Section IV. 



II. THE CHARM CONTENT OF THE POMERON 

After integration over the entire t range, the DDIS inclusive cross section can be written 

as 



d(3 dQ"^ dxw (3 



[l + {l-yf]F,^''>{(3,Q',x^), (1) 



where the contribution due to longitudinal structure function, F^^^\ has been neglected 
since it is expected to be small. Here the following kinematic variables are used to describe 
DDIS (in addition to usual DIS variables x, Q^, y, and W): 



and 



where Mx is the invariant mass of the diffractive system. The kinematical variable xp 
defined in Eq. (Q) can be interpreted as the fraction of the proton momentum transferred 
to the Pomeron, while /5, given by Eq. (^), may be considered as the momentum fraction 
of the Pomeron carried by the quark coupling to the photon. To simplify the notation, in 
what follows we will often write F2 instead of F2^*-^''(/3, Q^, a;p). 
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In a fit to the full data sample, HI Collaboration has found that a description of that 
considers only diffractive exchange requires a /^-dependent Pomeron intercept. However, this 
factorization breaking may be explained by introducing secondary trajectories [jl|. So, we 
suggest that, in the region where Pomeron exchange is the dominant process, the diffractive 
structure function could be expressed in a factorized form, 

F,'^{P,Q',xjp) = /i>/p(xp) F£{P,Q'), (4) 

where /f/p(xp) is the integrated Pomeron flux factor, and F^{i3, Q"^) is the Pomeron struc- 
ture function |2|,|l3. 

The contribution of 6-quarks to F2 is expected to be negligible due to the large mass of 
the bottom quark (as it takes place for DIS). Thus, we can omit this contribution and write 

F^{P,Q')=PY.^lF^\P^Q')^ (5) 

a 

Ca being the electric charge of the quark a (a = u, rf, s, c). 

For ^ ml, where ttIc is the charm-quark mass, we can regard u, d and s quarks to be 
massless and put (both the quark and antiquark are included in the distribution gp(/?, Q^)) 

Fi-\P,Q') = Fi'\p,Q') = F^\p,Q') = q^{P,Q'), (6) 

that results in 

F^{P,Q') = qAP,Q') + Fi^\p,Q',ml). (7) 

Recently, a factorization theorem has been proved for diffractive lepton scattering off 
nucleons [|I4[ from which structure functions of DDIS coincide with DIS structure functions. 
Therefore, quark and gluon distributions inside the Pomeron, gp(/?, Q^) and g^iP, Q"^), obey 
the same set of RG evolution equations as quark and gluon distributions inside the proton 



do. As the observed values of [3 are not too small, DGLAP equations can be used to 
perform such an evolution. 

In the present analysis, we suppose that charm quarks are mainly produced by virtual 
photon-gluon fusion and do not take part in the evolution of the light quarks. In such a 



case, by analogy with charm contribution to F2 0, we get the following equation for the 
charm contribution to DDIS structure function F2 , 

in which Qq = 2 GeV is assumed. 

Now, we isolate a similar term in q]p{f3, Q"^) and call the rest of Q"^) "direct contri- 
bution", that is 

,M Q'') - Q') + § 1; ^ c,(., Q^ P. 0) (f . P) . (9) 

Let us define the quantity 

AF^\^,Q^ml)= /f ^/'^[C,(^,g^P,0)-C,(;.,Q^P,m^)] 

By using these definitions, from Eqs. (0)-(|l3) we obtain 



+ 



2 



It follows from Eq. ([TT|) that 



Fif)(/3,Q^m,2) = ^ 



(12) 



An analogous difference between DIS structure functions with and without charm con- 
tribution, AF2{x , Q'^ , ml) , was calculated in Ref. [Q, where it was shown that it scales at 
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high Q"^ . The generahzation for DDIS is straightforward and the result (up to corrections 
0(m2/Q2)) reads 



AF(f)=A4^)(/3,m^) 



I k'^ Jp z yk"^ ' J dlnk'^ 



9 f(3 ,2' 



(13) 



In as order the expression for AC is of the formQ 



AC{v,u) = — \P,g{u) In 



7r 



1 + 



ml — u) 



2uf 



V + u(l — u) 



(14) 



Formula (|T^) does not contradict the factorization theorem for DDIS |14]. Namely, if we 
put A;^ = in Cg{P, Q^, k"^, ml) as it is usually done, the main contribution in Eq. (|) is due 
to the region /c^ ~ Q^, and one gets 



^ dz ^ I m:: 
a 



(15) 



On the other hand, in the difference of the diffractive structure functions, AF^\(3,ml^^ 
Eq. (0), the leading contributions cancell out. The quantity AC has the asymptotic be- 
havior 



(16) 



and the main contribution to the integral in k in Eq. (|13D comes from the region k 



~ m. 



i- 



As can be seen from Eq. (|T2D, Fp-* is defined via AFp'*. In its turn, AFp"* is given by 
formula (|13]) which contains a derivative of g-p in \nk^. This is related to the fact that we 



started from an exact expression for the coefficient function Cg depending on both 



and k"^. Thus, one can expect that the charm contribution to the DDIS structure function 



^ There are two misprints in formula (40) of Ref. [^j that were corrected in Eq. ([l^): the expression 
Pqg{y) = 1/2[(1 — y)^ + y^] should be put within the curly brackets and the factor as/47r should 
be replaced by ag/'i'K. 
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should significantly be dependent on both the form and evolution of the gluon distribution 
inside the Pomeron. 

For numerical estimates we shall use the quark and gluon distribution functions which 
have been obtained in Ref. by fitting the data on ^2^^^'' from HI and ZEUS collabora- 
tions |TB,[T^. Since in this analysis it was assumed Nf = 3 {Nj being the number of flavors), 
we have to rewrite Eq. ([12D in terms of corresponding parton distributions, gp"* and Let 
us also define gp"* (g^^) to be a quark (gluon) distribution for the case Nf = 4. 

It is useful to introduce the quantities 



— Q-p - Qip 



.(3) 



(17) 



and 



(18) 



It should be noted that 



Aqjp 



p • 



(19) 



From Eq. (|T2|), one obtains 



+ 



Nf=3 



+^r^^/>.p.(.)fA.p(f,.^) 



(20) 



9JqI P tt 

where the subscript |Af^=3 means that the corresponding quantities in the right hand side 
(RHS) of Eq. (|20D should be calculated with the use of the distributions qjp and g]p . 

The next step is to estimate the quantity A(yfp which enters into Eq. (|20|). Due to our 
assumption (no charm in light quark evolution), gp"* and gp"* obey one and the same DGLAP 
evolution equation [|T5|, 



+ 



+ 



dk^ as /•! dz 



Ql IT 



(21) 



for n = 3,4. The factor 3 in front of the first integral in the RHS of Eq. ( pT]) is related to 
the number of light flavors. 

As for initial quark and gluon distributions inside the Pomeron, we have gp\l3,Ql) 7^ 
g^^\f3,Ql), while 

Aqjp{(3,Ql) = 0, (22) 

that is no intrinsic charm in the Pomeron. 

If we neglect the variation of as with the change of the flavor number from Nf = 3 to 
Nf = 4, from Eq. we get 

The QCD-evolution parameter 

where 12Tth = 33 — 2A^j, rises slowly in and is numerically small even at rather high 
values of Q^. Por instance, for Qq = 2 GeV and A = 0.2 GeV we find ^(lO^ GeV^) ~ 0.13, 
.^(10^ GeV^) ~ 0.24. In particular, it enables one to solve the DGLAP equations by using 
an expansion in the parameter ^ |18[. 

From all said above, we obtain (up to small corrections 0(^^)) 

Due to the fact that Cg{z, Q^, k"^, ml) has no large logarithms, at k"^ ~ (see an exphcit 
expression for Cg{z, Q^, k'^, ml) in Ref. [^]), we obtain from Eq. (|^) in the leading logarithmic 
approximation (LLA) the expression 

where Pcg{z,m1/Q'^) = —dCg{z,Q'^,k'^,ml)/dlnk'^\i.2^Q2 is the modified form of the Pqg 
splitting function for the charm quark |jl9| . 



Now let us define 



1 dz 



P<ig{z) Ag 



'/3 



p 



(27) 



and 



(28) 



and put Tc 



r. This means that we ignore additional subleading m^-dependent terms 



in Peg with respect to Pgg. The quantity r may in principle depend on both P and Q^. 
Numerical estimates have shown, however, that it is weakly dependent on variable Q^. 
From Eqs. (|l^), (p5|)-p8|), we calculate r and rewrite F^^ in the following form 

2 AB 



2 _2\ ^ 

3A + B'' 



where 



(29) 



Nf=3 



and 



B 



a. 



I" dz Pgg{z) ^ g^^^ 

Jp z \z 



(30) 



(31) 



Here F^ is the Pomeron structure function, while q^"^ and AFp"* are defined in Eqs. (§) and 
(II)- 



III. THE CHARM CONTRIBUTION TO THE DIFFRACTIVE STRUCTURE 

FUNCTION 

In this section, we present quantitative results obtained in the present analysis as well 
as some comparison with other models. 
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A. Results of the present analysis 

Our concern now is the calculation of the charm contribution to F^^^\i3, Q^, xjp) in two 
different approaches and for different shapes of the quark and gluon distributions inside the 
Pomeron. In one approach the standard flux factor is employed, whereas in the other the 
renormalized flux factor is used (for brevity, we will refer to these quantities hereafter as 
STD and REN flux factors, respectively). For the former, it was assumed the Donnachie- 



Landshoff expression [^ , 



/5TD(xp,t) = g[Fi(t)pxlr^-(*) (32) 



while the latter is determined from the procedure prescribed in pT|, that is 



where 



r / ,x fsTD{xjp,t) , . 

W(xp,t)= ^^^^^^^^^^ (33) 



^(^Pmin) = / dxjp fsTD{xw,t) dt. (34) 

Jxtp . Jt=—oo 



By introducing Eq. ( p2D into Eq. (0) and assuming an exponential approximation for the 
form factor, -Ff (t) — 6^"*^*^ one obtains 

iV(xp_) = K^— [E,(7 - 2e Inxp^^J - ^,(7 - 2e lnxp_)], (35) 

where Ei{x) is the exponential integral, K = 9/9Q/47r^ and 7 = hot/a'. The minimum 
value of xp is x-p^^^ = {nip + m.^y/s for soft diffractive dissociation and x-p^^^^ = Q'^/f3s for 



DDIS [|T 



The distributions of the quarks and gluons inside the Pomeron, gp"* and were ob- 
tained from HERA data |16|,|T^ in Ref. [|T2| (we refer the reader to this paper for details). 
The parametrizations for each flux factor are described below. No sum rules were imposed 
on them to perform the fitting. 

Fit 1: Parametrizations obtained in with STD flux in which both quark and gluon distri- 
butions have a hard shape at the initial scale of evolution: 
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3/3 g^'^(/3,Q^) = 2.55/3 (1-/3), 

/3^7t'^(/3,Qo) = 12.08 /3 (1 -/3). (36) 

Fit 2: Parametrizations obtained witli tlie STD flux; tlie initial distributions correspond 
to a super-hard profile imposed to gluons by a delta function while quarks were left free to 
change according to the data: 

3/3 ggH/3,Q^) = 1.51 P'-'' (l-/3)°■«^ 

f3gS\f3,Ql) = 2.06 6(1 -P). (37) 

Fit 3: Parametrizations obtained with the REN fiux factor and a initial combination of the 
type hard-hard: 

3/3 gg^(/3,Q^) = 5.02/3 (1-/3), 

/3^7g^(/3,Q^) = 0.98 /3 (1 -/3). (38) 

All these three combinations of fiux factors and parton distributions of the Pomeron 
were applied in the calculation of the charm contribution to DDIS structure function. This 
quantity is given by the formula 

F^'\p,Q',x^) = /p/p(xp) F^\p,Q'), (39) 

where /f/p(xp) stands for the integrated (over t) fiux factors mentioned above and the 
charm structure function of the Pomeron is defined by Eqs. (|29D -(^ID. 

The results of our calculations are shown in Figs. 1-3. The upper curve in each figure 
corresponds to the total diffractive structure function, xjpF^, while the lowest one describes 
its charm component, X]pF2'^\ The difference X]p{F.p — Flf*) is also shown. 

In these figures, the theoretical results are presented together with recent HI and ZEUS 
data on |1],|2| which were not used in the fitting procedure mentioned above. The idea 
is not providing a precise description for these data, but giving the reader a possibility to 
compare the net charm contribution to the precision of present-day data. 
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As one can see in Figs. 1-2, the charm contribution to the diffractive structure function 
obtained with the STD flux factor amounts to 30% - 40%, depending on the values of 
f3,x-p, and Q^. To compare, the non-diffractive structure function F2 contains between 10% 
(low Q"^) and 30% (high Q^) of charm at small x. From these figures, we see that the 
charm contribution to F<P grows with the decrease of xp and is a little bit larger for the 
hard gluon distribution (Fig. 1) than it is for the super-hard gluons (Fig. 2). However, for 
both parametrizations (Fig. 1 and Fig. 2) it is comparable with the experimental errors^ of 
the HI and ZEUS data and, consequently, can likely be measured in forthcoming HERA 
experiments on diffractive dissociation processes. 

On the other hand, for the renormalized flux factor the charm component is very small in 
the full range of P and presented in Fig. 3. The reason is that the initial gluon distribution 
for this case, Eq. ([38|), is much smaller as compared to the initial gluon distribution with 
the same form for the standard flux factor, Eq. (^). 

Another way of comparing these results is shown in Fig. 4 in terms of F^\i3, Q"^), which 
is calculated for the three combinations of Pomeron flux factors with the respective structure 
functions considered here. 

Two main features that characterize our predictions for Fp\f3,Q'^) are evident in this 
figure: (1) the shape of the (3 distributions are quite similar (they are moderately hard at the 
initial scale) and change similarly with evolution; (2) the amount of charm is different 
in each case with the proportions seen in the figure. 

B. Comparison with other models 

The diffractive production of the open charm in DIS has been studied in the framework 
of perturbative two-gluon exchange between the cc-pair and the proton in Refs. |2^. In 
Refs. , non-perturbative approaches were used to calculate cross sections and spectra for 

^Statistical and systematic errors have been added in quadrature. 
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charm quark pair production. 

One common aspect of some of these models (the first two of Refs. and the first one 
of Ref. which is in contrast with the resuhs of our analysis shown in Fig. 4 is that their 
predictions for the charm contribution practically do not change at low (3 with evolution. 
Another distinctive feature of these models in respect to ours is that the (3 distributions are 
generally peaked at some intermediate (3 value that becomes larger with increasing Q^. This 



last aspect is also observed in the analysis by Levin et al. p2[ , although in this case the low 
j3 behavior does not follow the others. 

Another general observation is that the obtained steep rise of the charm component 



towards small is in qualitative accordance with the results of Refs. [22 



In Fig. 5, we present a quantitative comparison of our results for the charm contribution 
to F2^^\x]p, (3, Q^) with of those obtained by hotter |2^ for two values and xp = 0.001. 
It is seen that, in terms of the amount of charm, hotter 's predictions are comparable only 
to our renormalized case (Fit 3), although in terms of shape these distributions are quite 
different. Let us note, however, that hotter 's model is not adequate to describe diffraction 
in the complete /3-range as was mentioned by the author . 

Now let us consider other models, reminding that our analysis was performed in the con- 
text of the Ingelman-Schlein model. Predictions for the charm contribution to the Pomeron 



structure function have been made by using the same scheme in Refs. [^. However, no 



estimates of the charm contribution to the diffractive structure function F2 have been 



Ji UliC Uliailii UUlitllUUtlUli UllC UllllClUljlVC OUIUUUUIC lUliUUlUll -I 

presented. 



In Fig. 6, we present a comparison of our predictions for Fj^\j3, Q"^) with those obtained 



by Haakman et al. p4[. We see that in their analysis the charm structure function is pretty 
soft even at low where our results are predominantely hard. In terms of amount of charm, 
their results are comparable only to our Fit 3. 
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IV. CONCLUDING REMARKS 



We have considered in this paper the charm content of the Pomeron and its effects on 
the structure function measured in diffractive deep inelastic scattering. 

In the present analysis, the formulas are derived in a way to define this contribution from 
the quark and gluon distributions inside the Pomeron obtained previously by fitting HERA 
data on diffractive deep inelastic scattering. Two parametrizations have been chosen for the 
standard Pomeron flux factor corresponding to the hard and superhard gluon components 
of the Pomeron, whereas for the renormalized flux factor, the hard parton parametrization 
has been analyzed. 

Numerical calculations show that the results depend crucially on the Pomeron flux fac- 
tor. In particular, the charm content of the Pomeron is expected to be very small for the 
renormalized flux factor. As for charm contribution corresponding to the standard flux fac- 
tor, the estimates obtained allow us to think that it could be extracted from diffractive deep 
inelastic process with open charm production, taking into account the planned upgrades of 



the HERA experiment |25 
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Figure Captions 

Fig. 1 - Theoretical estimations of diffractive DIS structure functions in comparison with 
HERA data. The curves correspond to the total diffractive structure function, x-pF^ (solid 
curves), its charm component, a;pF2^'* (dotted curves), and the difference xtp{F^ — ^2"^'') 
(dashed curves). In the theoretical calculations were employed the flux factor and initial 
parton distributions of the Pomeron described as Fit 1 (see text). The experimental data 
are from HI IQ (filled circles) and ZEUS (open circles) collaborations. 

Fig. 2 - The same as Fig. 1, but with theoretical curves calculated from the formulas of Fit 

2 (see text). 

Fig. 3 - The same as Fig. 1, but with theoretical curves calculated from the formulas of Fit 

3 (see text). 

Fig. 4 - Predictions for the charm structure function in diffractive DIS as obtained with the 
parametrizations of Fit 1 (a). Fit 2 (b), and Fit 3 (c) and their respective evolution. 

Fig. 5 - Comparison of the charm contribution to F^^^\x-p, (3,Q'^) obtained in the present 
analysis (Fits 1, 2, and 3) with the predictions by hotter p2[ for two values and for 
xp = 0.001. 

Fig. 6 - Comparison of the charm structure function obtained in the present analysis (Fits 



1, 2, and 3) with the predictions by Haakman et al. p4 |. 
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Fig. 3: R.J.M.Covolan, A.V.Kisselev, M.S. Scares 
Title: Charm contribution to the Structure... 
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Fig. 5: R.J.M.Covolan, A.V.Kisselev, M.S.Soares 
Title: Charm contribution to the structure... 
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Fig. 6: R.J.M.Covolan, A.V.Kisselev, M.S.Soares 
Title: Charm contrilsution to the structure... 
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